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Abstract 

I review the current status of lattice Monte-Carlo computations of the 
equation of state of QCD, our current understanding of the thermaliza- 
tion process at collider energies, and two new theoretical developments, 
one pertaining to the ab-initio calculation of the low- a; gluon structure of 
heavy nuclei, the other to the chaotic behavior of QCD. At the end I give 
a brief overview of the status of various quark-gluon plasma signatures. 

1 Introduction 

In the preceding talk, Johanna StachelO surveyed the status of relativis- 
tic heavy ion experiments with fixed targets at the Brookhaven AGS and 
the CERN-SPS. Whereas the observed hadron yields and spectra pro- 
vide strong evidence that the final state — at the moment of break-up — is 
approximately thermal, we have little information about when the ther- 
malization is achieved during the course of the nuclear reaction. Since 
the energy deposition processes at AGS and SPS energies are dominated 
by soft strong interactions, we currently also lack the theoretical tools 
to make reliable calculations. Hence we do not know yet, and may never 
know, whether a thermalized quark-gluon plasma is formed in these re- 
actions. 

The situation appears to be much more favorable at the higher ener- 
gies where the future heavy ion colliders, RHIC and LHC, will operate. 
Two developments over the past few years have raised the hope that 
accurate calculations of the initial phase of nuclear reactions at such en- 
ergies will eventually be possible. Firstly, numerical simulations of lattice 
QCD are rapidly approaching the stage where quantitatively reliable pre- 
dictions of the equations of state of strongly interacting matter and its 
dynamical properties will be feasible. Secondly, the inclusion of color 
screening into parton cascades and nuclear structure functions opens up 
the prospect of a seamless description of the reaction dynamics within 
the framework of perturbative QCD from the initial nuclear structure 
up to a thermalized quark-gluon plasma. If this expectation comes to 
fruition, it will be possible to test quantitative theoretical predictions at 
RHIC and at the LHC. 

In this talk I will review the current status of lattice Monte-Carlo 
computations of the equation of state of QCD, our current understand- 
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ing of the thermalization process at collider energies, and two recent new 
theoretical developments, one pertaining to the ab-initio calculation of 
the low- a; gluon structure of heavy nuclei, the other to the chaotic be- 
havior of QCD. At the end I will give a brief overview of the status of 
various quark-gluon plasma signatures. 

2 The QCD Phase Diagram 

Rigorous results about the phase diagram of QCD are based on Monte- 
Carlo simulations of lattice-QCD at finite temperature. For the pure 
SU(3)-gauge theory, quantitatively reliable results including an extrap- 
olation to the continuum limit, have recently become available. Calcu- 
lations by the Bielefeld group □ on lattices up to 32'^ x 8 yield a critical 
temperature (in units of the string tension a) of 

Tc/aA = 0.629 ± 0.003, (1) 

or Tc = 260 MeV if the value of a derived from the Regge slope is 
used. The phase transition is of first order with a latent heat of slightly 
less than 2T^. The simulations provide clear evidence of significant in- 
teractions extending far above Tc- Fitting to predictions of thermal 
perturbation theory one finds, e.g., 

g^{5n) ^ 1.5. (2) 

The simulations also provide information about quantities of dynamical 
relevance, such as the speed of sound, which is found rising from < 0.1 
near Tc to cl ^ 1/3 at T = 5Tc. 

Simulations of full lattice-QCD with light dynamical quarks have not 
yet reached this advanced stage, but have also made significant progress. 
State-of-the-art calculations □ show a steep jump in the quantity e/T* 
at Tc ~ 150 — 200 MeV, indicating a rapid unthawing of new degrees 
of freedom. The uncertainty in the value of Tc is mainly caused by dis- 
crepancies between different ways of setting the scale (string tension or 
hadron masses). Although it remains unclear whether thermodynamic 
quantities will show a singularity at Tc in the continuum limit and what 
order phase transition it may be, the simulations clearly demonstrate 
that the quark condensate (i/'t/') drops steeply at exactly the same tem- 
perature: color deconfinement and chiral symmetry restoration occur at 
the same temperature. The results are expected to reach a similar qual- 
ity as those presently existing for the pure gauge theory when the next 
generation of parallel computers with close to teraflops performance will 
become available. 

Perturbative approaches, valid far above Tc, havsiiftlso made signif- 
icant progress. The equation of state is now known uij up to order g^, 
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and a scheme yielding a complete racrturbative result with only minimal 
lattice input has been worked outlj. The "bad news" is that it is now 
understood that the convergence of the perturbative series (even in the 
sense of an asymptotic series) requires g < 1 or as{T) < 0.1, which only 
holds for temperatures above the electroweak unification scale Q. 

Nevertheless, thermal perturbation theory provides important insight 
into the dynamics of the high-temperature phase of QCD: 

1. quarks and gluons develop dynamical masses of order gT. 

2. Long-range color forces (except static color magnetism) are dynam- 
ically screened on distances of order (gT)^^ , providing an effective 
infrared cut-off for many transport coefficients. 

3. Colored particles propagating through the plasma constantly ex- 
change their color with the medium, rendering the quark-gluon 
plasma a poor color conductor. As will be discussed below, this 
"color chaos" facilitates the rapid equilibration of degrees of free- 
dom in the plasma. 

3 Initial conditions at RHIC and LHC 

Most recent theoretical predictions for the initial conditions at which a 
thermalized quark-gluon plasma will be produced at heavy ion colliders 
are based on tke concept of perturbative partonic cascades. The parton 
cascade model □ starts from a relativistic transport equation of the form 

p^^F,{x,p)=a{x,p\Fk) (3) 

where Fi{x,p) denote the phase space distributions of QCD quanta. The 
collision terms d are obtained in the framework of perturbative QCD 
from elementary 2^2 scattering amplitudes allowing for additional 
initial- and final state radiation due to scale evolution of the perturba- 
tive quanta. To regulate infrared divergences, the parton cascade model 
requires a cut-off for the 2 — > 2 scattering (usually = 1.5 — 2 GeV/c) 
and a cut-off for time-like branchings (/^o — 0.5 — 1 GeV'^/c^). 

Numerical simulations of such cascades for heavy nuclei provide a 
scenario where a dense plasma of (predominantly) gluons develops in 
the central rapidity region betiween the two colliding nuclei shortly after 
the impact. Detailed studies u show that the momentum spectrum of 
partons becomes isotropic and exponential, i.e. practically thermal, at 
a time r « O.TAz in the rest frame of a slab of width at central 
rapidity. To allow for a hydrodynamic description, the width of the slab 
should exceed the mean free path of a parton. Including color screening 
effects, one finds that the mean free path of a gluon in a thermalized 
plasma is A/ « (3asT)~^ where T is the thermal slope of the parton 
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spectrum. For the very high initial values (T > 0.7 GeV) one concludes 
that a thermal hydrodynamic picture makes sense after Ti ~ 0.3 fm/c. 

The high density of scattered partons m A + A collisions makes it 
possible to replace the arbitrary infrared cut-off parameters p™'" and 
/io by dynamically calculated medium- induced cut-offslij. The dynam- 
ical density-dependent screening of color forces eliminates the need for 
p™'", and the suppression of radiative processes provided by the Landau- 
Pomeranchuk-Migdal effect makes the virtuality cut-off /xq superfluous. 
Note that the viability of this concept crucially depends on the high par- 
ton density: the dynamical cut-off parameters must lie in the range of 
applicability of perturbative QCD. Since the density of initially scattered 
partons grows as (Ai A2)^^'^(ln s)'^, this condition requires both large nu- 
clei and high collision energy. The calculations indicate that this crite- 
rion will be met at RHIC and LHC but not at the presently accessible 
energies of the SPS and AGS. The framework is also not applicable to pp 
or pp collisions at current energies because the parton density remains 
too low. 

Be dynamic screening of parton cascades can be implemented as fol- 
. According to the uncertainty principle, a parton-parton collision 
can be considered as complete after a "formation time" r/ ~ Ti/pr in the 
cm. frame, where pr is the momentum transfer. Accordingly, harder 
collisions are completed first. One can then consider the hard collisions 
as effectively screening the softer ones. The screening mass applicable 
to parton collisions of scale pr is then determined as 

„ POO 

l^'iPT) = -^Mp't) / d\\V, f{k)\ (4) 

where f{k) is the momentum density of more violently scattered partons. 
Feeding y?{pT) back into the differential cross sections determining the 
number of scatterings one obtains a differential equation for /(fc) or, 
equivalently, //^(pt). For Au + An collisions at RHIC energy, the screen- 
ing mass saturates at slightly below 1 GeV at small pr, and at 1.5 GeV 
for Pb -I- Pb collisions at the LHC. Both these values are comfortably 
within the range of applicability of perturbative QCD, demonstrating 
that there may be no need for an artificial infrared cut-off. The screen- 
ing of parton scattering by already scattered partons is analogous to the 
interactiori-among ladders in the traditional picture of soft hadronic in- 
teractionstJ. It would be interesting to rederive these results from this 
alternative point of view. 

The self-consistent parton cascade makes parameter-free (though 
somewhat model dependent) predictions about the initial conditions 
achieved at RHIC and LHC. For the heaviest nuclei one expects thermal- 
ization to occur at T = 730 MeV (RHIC) or Tc = 1150 MeV (LHC) and 
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initial energy densities of about 60 GeV/fm^ (RHIC) or 430 GeV/fm^ 
(LHC). At the initial moment (r ~ 0.25 fm/c) the parton plasma is not 
chemically equilibrated. It takes another several fm/c to achieve chemi- 
cal equilibrium. Following the evolution with the hydrodynamical model 
including self-consistent screening, the quark-gluon plasma is expected 
to last for about 5 fm/c until the temperature falls to Tc at RHIC and 
somewhat longer at the LHC. Final multiplicities are predicted to be 
dN/dy ^ 1700 (RHIC) and « 8000 (LHC). 

4 New Theoretical Approaches 

4-.1 Semiclassical parton structure 

While representing a major advance in the application of perturbative 
QCD to relativisitic nuclear collisions, the self-screened parton cascade 
still relies on the input of experimentally measured parton structure func- 
tions rfithe colliding nuclei. A new approach, due to the Minneapolis 
group t3 promises to make these structure functions themselves calcula- 
ble within perturbative QCD. The basic idea underlying this approach 
is that, as seen by partons with x < 10~^, the valence quarks in a heavy 
nucleus constitute a very dense, sheet-like random color source. The 
large area density p « 3A/nR^ defines a large scale parameter p? = p 
at which the QCD coupling ois{p^) is weak. Introducing light-cone coor- 
dinates it is then possible to formulate a systematic program for calcu- 
lating gluon and quark (sea) structure functions at small x as generated 
by classical random color fields and their quantum fiuctuations. At the 
classical level, corresponding to the Weizsacker-Williams approximation, 
the gluon density for not too small pr is given by 

1 dN ^ q.m'(A^c - 1) 

■kR^ dxdp"^ TT^XPy 

Loop corrections will introduce (lna;)/a; contributions, possibly leading 
to a power-like behavior at small x upon resummation. 

The approach can be extended to the collision between two nuclei, 
viewed as thp-interaction among two counter-propagating sheets of va- 
lence quarks tJ. Since screening effects are already included in this ap- 
proach, it may provide another description of the initial state produced 
in a nuclear collision at RHIC or LHC. 
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5 Chaotic Glue and Thermalization 

A completely different approach to the thermalization problem in non- 
Abelian gauge theories emerges from solutions of the clasical Yang-Mills 
equations. Numerical calculations have shown that these amisktions form 
a strongly chaotic, infinite-dimensional dynamical system lj'lj. The ev- 
idence is derived from the study of the (properly defined) distance be- 
tween two initially neighboring field configurations as a function of time. 
Neighboring configurations diverge at an exponential rate: 

Pi'W-A^WII-e^*. (6) 

The complete spectrum of Lyapunov exponents \i can be obtained from 
numerical solutions of the Yang- Mills equations on a lattice in Minkowski 
space. There is evidence that the lattice fields thermalize (in the micro- 
canonical sense) and the Lyapunov exponents scale as (g^T). 

A fundamental result of nonlinear dynamics is that the sum of all 
positive Lyapunov exponents (the so-called Kolmogorov- Sinai entropy) 
measures the rate of entropy growth after coarse graining. The Lya- 
punov exponents, therefore, provide a measure of the thermal equilibra- 
tion tii&fi. Using the largest Lyapunov exponent for the SU(3) gauge 
theory Lj, one finds a characteristic time of order 0.2 - 0.3 fm/c in the 
range of temperatures (300-700 MeV) relevant for RHIC. This agrees 
nicely with the results from parton cascade simulations. 

6 Quark-Gluon Plasma Signatures 

A wide variety of signatures for the formation of a quark-gluon plasma 
in heavy-ion collisions has been proposed. Considerable progress has 
been made iiv^ec^nt years in understanding the background to many of 
these signals lj'Ej. For example, it will be very difficult, if not impos- 
sible, to identify direct lepton pairs emitted by the plasmaj-because of 
a formidable background of decay leptons from D-mesons c3. On the 
positive side, the suppression of charmonium states observed in p + A 
collisions is becoming undsfstood as due to the absorption of the color 
octet component in nuclei cil. This provides a well-defined benchmark 
against which additional suppression effects can be identified. Indeed, 
the very high initial temperature now predicted for RHIC and the LHC 
will allow even charmed quarks to get into chemical equilibrium before 
hadronization. There is also a chance that even T mesons will be signifi- 
cantly suppressed at the LHC. The predicted production of multistrange 
baryons above thermal yields also remains a promising signature. 

It has recently been understood that the energy loss of ha*d partons 
in a quark-gluon plasma may be much higher than expectedEa. It may 
even grow with the length of the traversed medium. Accordingly, the 
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quenching of jet production, or of leading high-pT hadrons, has attracted 
increased interest as a probe of the early stage of the evolution of the 
plasma. 

The formation of chirally disorientated domains of the quark con- 
densate inpthe vacuum is a new, intriguing signature of the chiral phase 
transitioned. Such domains correspond to coherent excitations of the 
pion field; they would decay by producing large fluctuations away from 
1/3 in the tt^/tt ratio. Detailed numerical studies Ej of the linear sigma 
model have shown how disorientated domains can grow if the chiral tran- 
sition occurs rapidly, producing a temporarily unstable state of the chiral 
order parameter {if>'4^). Sophisticated techniemes for the identification of 
domain structures, such as wavelet analysisEj, may help observing these 
domains if they are produced experimentally. 



7 Summary 

As RHIC proceeds toward completion, theoretical advances during the 
last few years are giving us a much clearer view of the physics to be ex- 
pected in heavy ion collisions at that energy. Lattice-QCD simulations 
have unambiguously established the rapid "unthawing" of color and the 
restoration of chiral symmetry at a temperature below 200 MeV. Calcu- 
lations of parton transport processes based on perturbative QCD predict 
very high initial temperatures, above 500 MeV at RHIC and 1 GeV at 
the LHC. New approaches to low- a; parton structure and the thermal- 
ization problem hold the promise of a model and parameter-independent 
description of the formation and evolution of a quark-gluon plasma at 
RHIC energy and beyond. Finally, the various plasma signatures are 
becoming much better understood, providing valuable guidance for the 
experimental program at RHIC which will start in 1999. 
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